H=[1 -1;-1 2];
f=[-1;-4];
A=[1 1;-1 2;2 1];
b=[2;2;3];
lb=zeros(2,1);
 [x,fval,exitflag,output,lambda]=quadprog(H,f,A,b,[],[],lb)
